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Abstract 

Roughly one in three individuals is highly susceptible to motion sickness and yet the underlying causes 
of this condition are not well understood. Despite high heritability, no associated genetic factors have been 
discovered to date. Here, we conducted the first genome- wide association study on motion sickness in 80,494 
individuals from the 23andMe database who were surveyed about car sickness. Thirty- five single- nucleotide 
polymorphisms (SNPs) were associated with motion sickness at a genome- wide- significant level (p < 5 x 
10~^). Many of these SNPs are near genes involved in balance, and eye, ear, and cranial development (e.g., 
PVRL3, TSHZl, MUTED, H0XB3, H0XD3). Other SNPs may affect motion sickness through nearby 
genes with roles in the nervous system, glucose homeostasis, or hypoxia. We show that several of these 
SNPs display sex-specific effects, with as much as three times stronger effects in women. We searched for 
comorbid phenotypes with motion sickness, confirming associations with known comorbidities including 
migraines, postoperative nausea and vomiting (PONV), vertigo, and morning sickness, and observing new 
associations with altitude sickness and many gastrointestinal conditions. We also show that two of these 
related phenotypes (PONV and migraines) share underlying genetic factors with motion sickness. These 
results point to the importance of the nervous system in motion sickness and suggest a role for glucose 
levels in motion-induced nausea and vomiting, a finding that may provide insight into other nausea-related 
phenotypes such as PONV. They also highlight personal characteristics (e.g., being a poor sleeper) that 
correlate with motion sickness, findings that could help identify risk factors or treatments. 
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Introduction 

Motion sickness is provoked by exposure to a va- 
riety of motions (e.g., traveling in cars, boats, or 
planes; amusement park rides; skiing; and riding on 
camels) [1]. Simulators and virtual reality environ- 
ments can also induce motion sickness [2] . Symptoms 
of motion sickness include dizziness, nausea, vomit- 
ing, headache, and pallor [3]. Sweating, drowsiness, 
increased salivation, hyperventilation, and emotional 
distress may also occur. Motion sickness is associated 
with other conditions including migraines, vertigo, 
postoperative nausea and vomiting (PONV), and 
chemotherapy-induced nausea and vomiting (CINV) 
[1,4]. 

Roughly one in three individuals is highly suscepti- 
ble to motion sickness and the rest of the population 
may experience motion sickness under extreme condi- 
tions [5]. The underlying etiology of motion sickness, 
however, is not well understood. One theory suggests 
that motion sickness results from contradictory infor- 
mation the brain receives during motion [1,5]. The 
vestibular system of the inner ear, which senses mo- 
tion and body position and influences balance, signals 
"moving" to the brain, while the eye signals "sta- 
tionary" because the car or boat appears stationary 
relative to the viewer. The vestibular system is also 
thought to serve as a sensor of disequilibrium-causing 
neurotoxins (i.e., a toxin detector) and is believed to 
trigger the emetic response in order to rid the body 
of toxins. Thus, motion sickness may be an aber- 
rant trigger of the emetic response. Evidence for the 
involvement of the vestibular system comes from the 
observation that individuals with complete loss of the 
vestibular apparatus, a component of the vestibular 
system, are immune to motion sickness [1]. 

A variety of factors influence risk for motion sick- 
ness. Women are more susceptible than men [6-9] 
and younger individuals are at increased risk [8,9]. 
Ancestry may also play a role; there is some evi- 
dence that motion sickness occurs more frequently in 
individuals with Asian ancestry compared to Euro- 
pean ancestry [10,11]. Some variables are situational 
and/or behavioral. For instance, one study showed 
that passengers without a view of the road ahead 
were about three times more likely to experience ill- 



ness [8] and another report suggested that adopting a 
wider stance may reduce motion sickness [12]. There 
is also evidence that diet and eating behavior influ- 
ence risk [7]. 

Perhaps the most important and least understood 
variable is the underlying physiological susceptibil- 
ity of the individual. In women, increased Corti- 
sol levels are predictive of motion sickness [13] and 
susceptibility to motion sickness changes as a func- 
tion of the menstrual cycle, suggesting that lev- 
els of estrogen and other hormones might play a 
role [14]. In both sexes, hyperglycemia is impli- 
cated in motion-induced nausea and vomiting [15]. 
There is also some evidence that lower baseline lev- 
els of adrenocorticotropic hormone (ACTH) [16], also 
known as corticotropin, and low sympathetic nervous 
system activity [17] increases susceptibility. Finally, 
since antihistamines (e.g., Dramamine), anticholiner- 
gics (e.g., scopolamine), and sympathomimetics (e.g., 
d-amphetamine and ephedrine) are effective treat- 
ments, altered baseline activity of the receptors these 
drugs bind to might influence risk for motion sickness. 

Although heritability estimates for motion sick- 
ness range from 57-70% [18], genome- wide association 
studies (GWAS) on this phenotype have not been re- 
ported. Here, we describe a large GWAS in which we 
find 35 regions significantly associated with motion 
sickness. 

Results 

Genome- wide Association Study of Mo- 
tion Sickness 

We performed a GWAS in 80,494 individuals from the 
customer base of 23andMe, Inc., a personal genetics 
company. Participants were of primarily European 
ancestry and were at most distantly related to each 
other (i.e., first cousins and closer were excluded). 
Motion sickness was assessed using online self-report. 
Participants responded to questions about their de- 
gree of car sickness and questions were combined onto 
a scale of 0 (never motion sick), 1 (occasionally), 2 
(sometimes), or 3 (frequently). Details about the co- 
hort can be found in Table 1 and in the Methods. All 
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analyses were controlled for age, sex, and five princi- 
pal components of genetic ancestry. Manhattan and 
quantile-quantile plots are provided in Figures 1 and 
SI. 

Lead SNPs with p-values under 5 x 10~^ for mo- 
tion sickness are shown in Table 2; 35 regions were 
significant (Figure ??). We created a genetic propen- 
sity score based on the number of risk alleles for the 
35 index SNPs. Individuals in the top five percent of 
the distribution (allele dosage of 40.25 or more risk 
alleles) had an average motion-sickness score 0.546 
units higher than those in the bottom five percent 
(28.37 or fewer risk alleles). The top five percent 
had 6.37 times increased odds of being "frequently" 
motion sick as opposed to "never" motion sick as 
compared to the bottom five percent. The variance 
in motion sickness explained by the propensity score 
(which may be infiated as it was assessed in the dis- 
covery population) was 0.029. 

A few associated SNPs are in regions implicated in 
eye and ear development or balance. For example, 
our most significant association is with rs66800491 
{p = 4.2 X 10""^^), located roughly one Mb upstream 
of PVRL3, which encodes the cell adhesion protein 
Nectin-3. Loss of PVRL3 expression in both hu- 
mans and mice results in ocular defects [19]. The 
SNP rsl0514168 [p = 2.7 x 10"^) is located down- 
stream of TSHZl, a gene involved in inner ear de- 
velopment in the mouse [20]. Another association is 
with rs2153535 {p = 2.7x 10~^^), located upstream of 
MUTED^ which is implicated in balance [21]. Three 
additional associated SNPs are near genes with major 
roles in early development: rs2551802 (p = 2 x 10~^^) 
between H0XD3 and H0XD4 ; rs9906289 {p = 6.4 x 
10-^^) in H0XB3; and rsl49951341 {p = 3.4 xlO-^^) 
near TLEl The HOXD SNP is in LD (r^ ^ 0.9) 
with rs2072590, which is associated with ovarian can- 
cer [22]. 

Several other associated SNPs are located near 
genes involved in neurological processes includ- 
ing synapse development and function: rsl 1713169 
{p = 5.9 X 10~^^) in NLGNl encoding neuroligin; 
rs6069325 {p = 7.2 x 10"^^) upstream of CBLN4 en- 
coding a member of the cerebellin precursor protein 
family; rs62018380 {p = 1.7 x 10"^) downstream of 
MCTP2^ a gene involved in intercellular signal trans- 



duction and synapse function; rs7957589 in PDZRN4 
{p = 7.9 X 10-^^) near CNTNl (contactin 1), which 
plays a role in axon guidance during neural devel- 
opment [23]; and two independent SNPs, rs4343996 
and rs34912216 {p = 8.7 x 10-^^ and 2.7 x 10-^ 
respectively) in SDKl encoding sidekick- 1, a cell ad- 
hesion molecule that localizes to synapses. The SNP 
rs2150864 {p = 6.3 x 10"^^) is located about 1.5 Mb 
upstream of LING 02, a gene implicated in essential 
tremor [24]. Additional associated SNPs in or near 
genes in neurological pathways include: rs9834560 
{p = 9.7 X 10"^^) in CPNE4 encoding copine-4, and 
two independent SNPS in or near AUTS2 (rsll95218 
and rs6946969 {p = 4.5 x lO"^^ and 1.9 x 10"^, re- 
spectively) . 

Other associated SNPs are in regions involved in 
glucose and insulin homeostasis. For example, the 
second most significant association we found is with 
rs56051278 {p = 1.5 x 10-^9) in GPD2 that en- 
codes glycerol- 3-phosphate dehydrogenase 2, an en- 
zyme implicated in glucose homeostasis. This SNP 
is in high LD (r2 ^ 0.8) with rs2116665 (the non- 
synonymous substitution H264R in GPD2) that was 
previously associated with free fatty acid and glycerol 
levels [25]. The SNP rslll29078 {p = 3.4 x 10"^^) 
is located downstream of UBE2E2, which encodes a 
component of the ubiquitin-proteasome system. This 
system is implicated in the autophagy of pancre- 
atic beta-cells that produce insulin and plays impor- 
tant roles in insulin homeostasis [26]. In addition, 
rs705145 {p = 1.4 x 10~^^) is located just upstream 
of GPR26, encoding a G protein-coupled receptor. 
Mice with a deletion of the GPR26 gene develop 
hyperphagia and diet-induced obesity, which leads 
to metabolic complications linked to obesity includ- 
ing glucose intolerance, hyperinsulemia and dyslipi- 
demia [27]. The SNP rs4076764 [p = 2.9 x 10"^) 
is located upstream of RGS5, a regulator of G pro- 
tein signaling. Loss of RGS5 in the mouse is also 
associated with hyperphagia [28]. Finally, rs71 70668 
(p = 1 X 10"^^) is located upstream of NR2F2 encod- 
ing COUP-TFII (chicken ovalbumin upstream pro- 
moter transcription factor II), a protein with roles in 
glucose homeostasis and energy metabolism [29]. 

The remaining associated SNPs are in regions im- 
plicated in hypoxia (rsl858111 near RWDD3, p = 
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Figure 1: Manhattan plot. The 35 genome- wide significant regions are listed with the proposed candidate 
gene; regions that are close together share a label. 



Group 


Total 


Male 


Female 


< 30 


31-45 


46-60 


> 61 


Never 


40,042 


25,137 


14,905 


5,510 


10,707 


10,592 


13,233 


Occasionally 


24,902 


11,855 


13,047 


4,597 


8,451 


6,459 


5,395 


Sometimes 


6,723 


3,067 


3,656 


1,175 


2,015 


1,759 


1,774 


Frequently 


8,827 


3,011 


5,816 


1,784 


3,173 


2,310 


1,560 


Total 


80,494 


43,070 


37,424 


13,066 


24,346 


21,120 


21,962 



Table 1: Cohort statistics for motion sickness GWAS. Degree of motion sickness stratified by sex and 
age. Females and younger people tend to be more motion sick. 



4.1 X 10 ^^); iron homeostasis (rsl0970305 near 
ACOl^ p = 1 X 10~^^); brown adipose tissue 
(rs61759167 in PRDM16, p = 3.5 x 10"^^); and 
other less characterized processes: rs2360806 {p = 
7.2x10-1^) in ST18, rs2318131 {p = 3.3x10-^2) near 
C0PS8, rs60464047 {p = 6.7 x 10"^) in POU6F2, 
rsl0752212 {p = 1.1 x lO'^^) near CELF2, rs6833641 
{p = 1.8 X 10-^) near ARAP2, rsl7515225 and 
rs34311235 (independently associated, p = 2.5 x 10~^ 
and 7.9 x 10"^, respectively) in LRPIB, rsl378552 
{p = 4.3 X 10~^) in a gene desert on 4q34.3; and 
rsl782032 {p = 9x 10"^) near TUSCl, and rsl847202 
near SHQl and GXYLT2. Finally, rs997295 in 
MAP2K5 {p = 3.3 x 10"^) is in LD with rs2241423 



(r^ ^ 0.36), which is associated with BMI [30]. 



Enrichment 

Analysis of all regions with p < 10~^ using GREAT 
[31] showed a significant enrichment in regions con- 
taining genes involved in fusion of atlas and occip- 
ital bones (FDR=0.002) and abnormal arcus ante- 
rior morphology (FDR=0.038) in mouse. The genes 
annotated with one or both of these processes were 
HOXB, HOXD, TSHZl, and RARE regions (the SNP 
near RARE is rs2067120, p = 8.2 x 10"^). 
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SNP 


band 


position 


alleles 


p- value 


effect 


95% CI 


Frequency 


Quality 


Gene 


rs66800491 


3ql3.13 


109,634,127 


A/G 


4.2 X 


10~ 


44 


-0.078 


(-0.089,-0.067) 


0.683 


0.982 


PVRL3 


rs56051278 


2q24.1 


157,381,754 


A/G 


1.5 X 


10~ 


29 


0.066 


(0.055,0.078) 


0.265 


0.992 


GPD2 


rsl0970305 


9p21.1 


31,372,583 


A/C 


1.0 X 


10~ 


27 


-0.057 


(-0.068,-0.047) 


0.505 


0.956 


ACOl 


rsll95218 


7qll.22 


68,624,342 


A/G 


4.5 X 


10~ 


22 


0.095 


(0.076,0.114) 


0.916 


1.000* 


AUTS2 


rs705145 


10q26.13 


125,226,178 


A/C 


1.4 X 


10~ 


21 


-0.051 


(-0.062,-0.041) 


0.638 


0.997 


GPR26 


rslll29078 


3p24.3 


22,592,321 


A/G 


3.4 X 


10~ 


21 


0.057 


(0.045,0.069) 


0.751 


0.984 


UBE2E2 


rs6069325 


20ql3.2 


54,139,486 


G/T 


7.2 X 


10~ 


21 


0.069 


(0.054,0.083) 


0.841 


0.933 


CBLN4 


rs2153535 


6p24.3 


8,369,679 


C/G 


2.7 X 


10~ 


18 


0.046 


(0.035,0.056) 


0.480 


0.969 


MUTED 


rs2150864 


9p21.1 


29,363,265 


A/G 


6.3 X 


10~ 


15 


0.042 


(0.032,0.053) 


0.348 


1.000 


LING02 


rs9834560 


3q22.1 


131,716,105 


A/C 


9.7 X 


10~ 


15 


-0.041 


(-0.051,-0.031) 


0.610 


1.000 


CPNE4 


rsl858111 


lp21.3 


96,089,731 


A/G 


4.1 X 


10~ 


14 


-0.039 


(-0.050,-0.029) 


0.567 


0.999 


RWDD3 


rs61759167 


lp36.32 


3,091,587 


C/T 


3.5 X 


10~ 


13 


0.047 


(0.034,0.059) 


0.231 


0.918 


PRDM16 


rsll713169 


3q26.31 


173,384,589 


A/C 


5.9 X 


10~ 


13 


-0.052 


(-0.067,-0.038) 


0.160 


0.918 


NLGNl 


rs2551802 


2q31.1 


177,022,158 


C/G 


2.0 X 


10~ 


12 


0.040 


(0.029,0.052) 


0.697 


0.951 


HOXD 


rs2318131 


2q37.3 


237,933,966 


A/C 


3.3 X 


10~ 


12 


0.038 


(0.027,0.049) 


0.343 


0.999 


COPS8 


rsl49951341 


9q21.31 


81,268,149 


A/C 


3.4 X 


10~ 


12 


-0.050 


(-0.064,-0.036) 


0.798 


0.794 


TLE4 


rs9906289 


17q21.32 


46,644,677 


C/T 


6.4 X 


10~ 


11 


0.083 


(0.058,0.108) 


0.046 


0.944 


HOXB 


rs2360806 


8qll.23 


53,125,734 


A/C 


7.2 X 


10~ 


11 


0.047 


(0.033,0.061) 


0.162 


0.958 


ST18 


rs4343996 


7p22.2 


3,362,642 


A/G 


8.7 X 


10~ 


11 


0.034 


(0.023,0.044) 


0.451 


0.994 


SDKl 


rs7170668 


15q26.2 


96,014,143 


C/T 


1.0 X 


10~ 


10 


0.035 


(0.024,0.045) 


0.632 


1.000 


NR2F2 


rsl0752212 


10pl4 


10,917,121 


A/G 


1.1 X 


10~ 


10 


0.034 


(0.024,0.044) 


0.532 


0.974 


CELF2 


rs7957589 


12ql2 


41,874,282 


A/T 


7.9 X 


10~ 


10 


-0.047 


(-0.061,-0.032) 


0.146 


0.930 


CNTNl 


rs62018380 


15q26.2 


95,275,917 


A/C 


1.7 X 


10' 


-9 


0.047 


(0.032,0.062) 


0.869 


0.960 


MCTP2 


rs6833641 


4pl5.1 


35,563,786 


C/G 


1.8 X 


10' 


-9 


0.046 


(0.031,0.062) 


0.852 


0.870 


ARAP2 


rs6946969 


7qll.22 


70,211,027 


A/G 


1.9 X 


10' 


-9 


0.033 


(0.022,0.043) 


0.658 


0.996 


AUTS2 


rsl7515225 


2q22.1 


141,545,755 


C/T 


2.5 X 


10' 


-9 


0.032 


(0.021,0.042) 


0.445 


0.958 


LRPIB 


rsl0514168 


18q22.3 


73,098,949 


A/C 


2.7 X 


10' 


-9 


-0.047 


(-0.062,-0.031) 


0.854 


0.861 


TSHZl 


rs4076764 


lq23.3 


163,441,286 


C/T 


2.9 X 


10' 


-9 


0.033 


(0.022,0.044) 


0.649 


0.933 


RGS5 


rs997295 


15q23 


68,016,343 


G/T 


3.3 X 


10' 


-9 


-0.033 


(-0.044,-0.022) 


0.588 


0.994* 


MAP2K5 


rsl378552 


4q34.3 


180,356,846 


C/T 


4.3 X 


10' 


-9 


-0.032 


(-0.043,-0.022) 


0.322 


0.998 


AG A 


rs60464047 


7pl4.1 


39,418,538 


A/T 


6.7 X 


10' 


-9 


-0.043 


(-0.057,-0.028) 


0.850 


0.956 


POU6F2 


rs34311235 


2q22.2 


142,767,433 


C/T 


7.9 X 


10' 


-9 


0.032 


(0.021,0.042) 


0.368 


0.955 


LRPIB 


rsl782032 


9p21.2 


25,804,285 


A/G 


9.0 X 


10' 


-9 


-0.031 


(-0.042,-0.020) 


0.544 


0.920 


TUSCl 


rsl847202 


3pl3 


72,934,371 


C/T 


2.5 X 


10' 


-8 


0.031 


(0.020,0.042) 


0.644 


0.945 


GXYLT2 


rs34912216 


7p22.2 


4,118,377 


A/G 


2.7 X 


10' 


-8 


-0.035 


(-0.047,-0.023) 


0.727 


0.849 


SDKl 



Table 2: Genome- wide significant index SNPs. Alleles are reported in alphabetical order with respect 
to the positive strand of build 37 of the human genome. The effect is the change per copy of the second 
allele on a four point scale of increasing motion sickness. Frequency is the frequency of the alphabetically 
second allele in the cohort. Quality is imputation for imputed SNPs, call rate for genotyped SNPs (those 
marked with a *). Gene is a proposed candidate gene in the region. 



Phenotypic Study of Motion Sickness 

We investigated comorbidities with motion sickness 
within the 23andMe database. Briefly, we looked 
at partial correlations between each of 695 different 
phenotypes and motion sickness, controlling for age, 
sex (where applicable), and 5 principal components. 
Table 3 shows selected large correlations. Some of 
the associated phenotypes are known symptoms of 
motion sickness (e.g., headache) or established co- 
morbidities (migraines, vertigo, PONV, and morning 



sickness). In addition to PONV, other gastrointesti- 
nal (GI) phenotypes were also associated with motion 
sickness (e.g., irritable bowel syndrome (IBS); acid 
reflux; stomach upset with antidepressants, codeine 
and nonsteroidal anti-inflammatory drugs (NSAIDs); 
and indigestion with dairy products). Other associ- 
ations include poor sleep, poor circulation, altitude 
sickness, hay fever, and neuroticism. Phenotypes as- 
sociated with lower risk for motion sickness include 
a history of tobacco use, a good sense of direction, 
higher BMI, being single, and a better ability to han- 
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die stress. 



Phenotype 


P 


AT 
iV 


Dizziness 


0.119 


24606 


PONV 


0.117 


25223 


Lightheaded during exercise 


0.114 


23434 


Vomiting from codeine 


0.114 


12176 


Altitude sickness 


0.111 


41661 


Morning sickness 


0.108 


13285 


Daytime sleepiness 


0.080 


29384 


Indigestion due to dairy products 


0.071 


29659 


Hay fever 


u.uoy 


OOO/I o 

zzz^z 


Headache after red wine 


U.UDD 




Vertigo 


u.uoo 


54934 


Back pain frequency 


0 065 


zooUo 


Neurotic ism 


0.063 


38711 


Rosacea 


U.UDo 




Irritable bowel syndrome (IBS) 


U.UDZ 




]VIosc[uito bites itching more 


0.061 


51292 


Greater perceived stress 


U.UDU 




IVIore colds last year 


0.060 


34759 


Drowsiness from Benadryl 


0.055 


19307 


Migraines 


0.055 


72901 


Seasonal affective disorder 


0.051 


32273 


GI pain from NSAIDs 


0.047 


29846 


More sleep needed 


0.046 


39115 


Nausea from antidepressants 


0.046 


12079 


Dry mouth from antihistamines 


0.033 


15826 


Acid reflux 


0.033 


68998 


Being single 


-0.033 


65186 


Adventurous 


-0.033 


34532 


BMI 


-0.034 


75217 


Punctuality 


-0.035 


37486 


Drive through food frequency 


-0.035 


41467 


Good sense of direction 


-0.039 


46440 


Positive attitude towards self 


-0.040 


34082 


Pack years (cigarettes) 


-0.046 


66381 


Sound sleeper 


-0.047 


48339 


Ability to handle stress 


-0.072 


30965 



Table 3: Selected partial correlations with mo- 
tion sickness. Partial correlations are controlled for 
age, sex, and 5 PCs. TV refers to the number of peo- 
ple with data for both motion sickness and the second 
trait. Traits are sorted by partial correlation. 



Genetic Correlations Between Motion 
Sickness and Related Phenotypes 

We determined if any of the 35 SNPs associated with 
motion sickness were also associated with six corre- 
lated and clinically important phenotypes (PONV, 
migraines, hay fever, altitude sickness, morning sick- 
ness, and vertigo). Table 4 shows SNPs associated 
with these phenotypes with a (Bonferroni-corrected) 
p- value under 0.05/35 ^ 0.0014 (under a more strin- 



gent threshold of 0.05/(35 x 6) ^ 0.0002 only the 
first two are significant). One motion sickness- 
associated SNP was significantly associated with mi- 
graines: rs61759167 in PRDM16 {p = Llx"^). A 
previous study [32] reported an association between 
migraines and another SNP in PRDM16^ rs2651899, 
which is in weak LD with rs61759167 (r^ ^ 0.44). 
Three motion sickness-associated SNPs were also sig- 
nificantly associated with PONV: rs6833641 near 
ARAP2, rsll95218 near AUTS2, and rs6069325 near 
CBLN4. For all four examples, the higher risk allele 
for migraines or PONV is also the higher risk allele 
for motion sickness. We did not detect significant as- 
sociations between motion sickness-associated SNPs 
and altitude sickness, hay fever, morning sickness, or 
vertigo. While this data suggests some shared etiol- 
ogy for motion sickness and PONV or migraines, it is 
difficult to assess whether or not this is due to shared 
causal SNPs. 

Sex-specific effects 

Motion sickness is much more common in women 
than in men (Table 1) and several of our SNPs show 
much stronger effects in women than in men. The 
SNP rs66800491 has a 1.5x larger effect in women (- 
0.097 versus -0.062) and rsl847202 has a 3x larger 
effect in women (0.048 versus 0.016) (both SNPs 
p < 0.05 for interaction, corrected for 35 tests). Over- 
all 26 of the 35 SNPs have estimated larger effects in 
women than men (binomial p < 0.003; Table SI). 

Discussion 

Here we report 35 novel, genome- wide significant as- 
sociations for motion sickness (Table 2). Genes in 
regions associated with motion sickness appear to 
play roles in eye and ear development, balance and 
other neurological processes, and glucose homeosta- 
sis. Two of the genome-wide significant regions con- 
tain hypoxia-inducible genes. We also provide evi- 
dence that motion sickness is phenotypically associ- 
ated with numerous conditions and traits (Table 3). 

Since motion sickness is thought to stem from the 
brain receiving contradictory signals from the inner 
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Phenotype 


SNP 


TV 


pvalue 


effect 


95% CI 


Migraine 


rs61759167 


72,901 


1.1 X 10-^ 


0.08 


(0.051, 0.119) 


PONV 


rsll95218 


25,223 


0.00012 


-0.14 


(-0.213, -0.069) 


PONV 


rs6069325 


25,223 


0.00079 


0.09 


(0.038, 0.143) 


PONV 


rs6833641 


25,223 


0.00101 


0.09 


(0.037, 0.148) 



Table 4: Significant associations between motion sickness-associated SNPs and other pheno- 
types. N is the number of people with data for motion sickness and migraines or PONV. 



ear versus the eye (e.g., the inner ear signals "mov- 
ing" while the eye signals "stationary"), it is inter- 
esting that a region implicated in eye development 
(rs56100358 near PVRL3) is our most significant as- 
sociation. Chromosomal rearrangements that lead to 
loss of PVRL3 expression have been associated with 
ocular defects in humans and the PVRL3 knockout 
mouse exhibits lens and other vision problems [19]. 
The associations with regions involved in the inner 
ear (rsl2111385 near MUTED and rsl435985 near 
TSHZl) are also interesting since disturbances in 
the vestibular system of the inner ear, which senses 
motion and body position and infiuences balance, 
are thought to play a central role in motion sick- 
ness. It has been suggested that the mouse homolog 
of MUTED controls the synthesis of otoliths of the 
vestibular labyrinth of the inner ear [21]. Otoliths are 
sensitive to gravity and linear acceleration and play 
a role in balance. Mutations in TSHZl and deletions 
in the 18q22.3 region that includes TSHZl are asso- 
ciated with congenital aural atresia (CAA) [33,34], 
a spectrum of ear deformities that involve malforma- 
tion of the external auditory canal (EAC). More gen- 
erally, our enrichment analysis suggests that genes 
involved in certain aspects of cranial developmental 
may play an important role in motion sickness. As- 
sociations with SNPs in or near genes involved in 
synapse formation and function {NLGNl, CBLN^-, 
MCTP2, PDZRNl CNTNl, and SDKl) and other 
neurological pathways {LING02, CPNE^ AUTS2) 
point to the importance of the brain in motion sick- 
ness. 

Five associated SNPs are in or near genes impli- 
cated in glucose and insulin homeostasis or BMI. 
Although these SNPs are not in linkage disequilib- 



rium (LD) with SNPs reported in GWAS of type 
2 diabetes (T2D) [35-40], rs56051278 is in high LD 
(r^ ^ 0.8) with rs2116665, a nonsynonymous sub- 
stitution (H264R) in the GPD2 gene. H264R has 
been associated with increased plasma glycerol and 
free fatty acid concentrations in a French Canadian 
population [25]. Increased free fatty acid levels are 
indicative of glucose intolerance and hyperinsuline- 
mia. Although it is unclear why genes involved in 
glucose and insulin regulation might also play a role 
in motion sickness, one study suggested that hyper- 
glycemia may be related to the gastrointestinal symp- 
toms of motion sickness [15]. In this study, individu- 
als who experienced motion-induced nausea and vom- 
iting had lower levels of insulin than people who did 
not experience gastrointestinal symptoms. The study 
further suggested that stable glucose levels might 
help to relieve motion-induced gastrointestinal upset. 

At least two of our associated SNPs are near 
hypoxia-inducible genes: RGS5 and RWDD3 (en- 
coding the RSUME protein). RSUME promotes the 
activity of hypoxia-inducible factor 1 (HIF-l alpha), a 
master regulator of the hypoxic response [41]; RGS5 
is an apoptotic stimulator induced by hypoxia in en- 
dothelial cells [42]. This data suggests a potential 
relationship between motion sickness and hypoxia. 
Motion sickness might lead to hypoxia or individuals 
predisposed to hypoxia might also be more suscepti- 
ble to motion sickness. Both possibilities are intrigu- 
ing since our phenotypic analysis suggested an as- 
sociation between motion sickness and altitude sick- 
ness, which occurs when individuals become hypoxic 
at higher altitudes (Table 3). 

Among the regional association plots (Supple- 
ment), one SNP in particular stands out: rsll95218 
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near AUTS2. This genotyped SNP has a value 
under 10~^^ and no other SNPs in the region have 
p < 10~^. This lack of signal from LD is not terribly 
surprising, as none of the three proxy SNPs (r^ > 0.2) 
for this SNP in 1000 Genomes pass imputation qual- 
ity control in our data. As the clusters for this SNP 
look excellent and the call rate is 99.98%, we believe 
this is a true signal. 

Certain phenotypic associations are interesting 
given what is known about motion sickness. PONV 
is an established comorbidity of motion sickness [43] 
and is thought to stem from the anaesthetics that are 
administered for surgery. It may not, therefore, be 
surprising that motion sickness was also associated 
with vomiting and/or nausea with use of codeine, 
antidepressants and nonsteroidal anti-inflammatories 
(NSAIDs). Additional gastrointestinal phenotypes 
(e.g., IBS, acid reflux, and indigestion with dairy 
products) as well as other drug-related phenotypes 
like being drowsy when taking Benadryl and feeling 
jittery when taking Sudafed were also associated with 
motion sickness. Interestingly, our findings suggest 
shared genetic susceptibility for both motion sickness 
and PONV (Table 4). 

Some phenotypic associations might provide clues 
about the etiology of motion sickness (e.g., poor cir- 
culation and becoming light headed with exercise) 
or they might suggest simple remedies for motion 
sickness such as improving sleep quality. A number 
of associated phenotypes were related to personality 
(e.g., neuroticism) or behavior (e.g., smoking). We 
note, however, that it is difficult to assess causal- 
ity for these phenotype-phenotype associations. For 
example, does being a sound sleeper make one less 
susceptible to being motion sick, or vice versa, or 
are both related to a third condition? The valid- 
ity of these novel phenotypic findings is bolstered 
by the fact that we also detected associations with 
known symptoms (dizziness and headache) and es- 
tablished comorbidities (PONV, migraines, vertigo 
and morning sickness) of motion sickness. In some 
cases we even identified shared genetic factors for mo- 
tion sickness and related comorbidities (e.g., PONV 
and migraines). Some of the correlated phenotypes 
are not established comorbidities or symptoms of mo- 
tion sickness, however, and do not have an obvious 



biological relationship to motion sickness. 

Our web-based method of capturing phenotypic in- 
formation allows us to build a very large cohort (e.g., 
80,494 individuals in our study), but we may not have 
obtained a complete picture of an individual's sus- 
ceptibility to motion sickness. For finding SNPs, the 
gain in power from having a large sample more than 
makes up for the reduction in power due to possi- 
ble misclassification. An additional potential limita- 
tion is that we only surveyed individuals about car 
sickness; future studies should investigate these SNPs 
in populations phenotyped for other forms of motion 
sickness. An advantage of our web-based phenotypic 
collection method is that we can easily investigate 
whether seemingly related traits have shared under- 
lying genetics. We identified four SNPs simultane- 
ously associated with motion sickness plus PONV or 
migraines. These findings may provide clues into the 
etiology of all three conditions and may point to over- 
lapping risk factors or treatments. 

Methods 
Human Subjects 

All participants were drawn from the customer base 
of 23andMe, Inc., a consumer genetics company. This 
cohort has been described in detail previously [44,45]. 
Participants provided informed consent and partici- 
pated in the research online, under a protocol ap- 
proved by the external AAHRPP-accredited IRB, 
Ethical & Independent Review Services (E&I Re- 
view). 

Phenotype collection 

Participants invited to fill out web-based question- 
naires, which included four questions about motion 
sickness during road travel, whenever they logged 
into their 23andMe accounts. The responses to each 
question were translated into a motion sickness score 
of 0 (never), 1 (occasionally), 2 (sometimes), or 3 
(frequently). Responses of 'Tm not sure" or "Don't 
know" were excluded from the analysis. The ques- 
tions were prioritized (1-4) and a participant's final 
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score was based on the question they answered that 
had the highest priority. The questions were: 

1. "How often do you experience motion sickness 
while in a car? (Never / Occasionahy / Some- 
times / Frequently / I'm not sure)" 

2. "Have you experienced motion sickness while 
riding in a car (car sickness)? (Yes, I do now 
frequently / Yes, I did frequently, but only as a 
child / Yes, occasionally / No / Don't know)" 

3. "As a child, how often did you experience motion 
sickness while in a car? (Never / Occasionally / 
Sometimes / Frequently / Fm not sure)" 

4. "Can you read in a moving car without becoming 
nauseated? (Never / Sometimes / Always / Fm 
not sure)" 

Responses to the final question were scored as 
Never=3, Sometimes=l, and Always=0. 

Genotyping and imputation 

Participants were genotyped and additional SNP 
genotypes were imputed against the August 2010 re- 
lease of the 1000 Genomes data as described pre- 
viously [46]. Briefiy, they were genotyped on at 
least one of three genotyping platforms, two based 
on the Illumina HumanHap550+ BeadChip, the third 
based on the Illumina Human OmniExpress+ Bead- 
Chip. The platforms included assays for 586,916, 
584,942, and 1,008,948 SNPs, respectively. Geno- 
types for a total of 11,914,767 SNPs were imputed 
in batches of roughly 10,000 individuals, grouped by 
genotyping platform. Imputation was performed as 
in [46]. Prior to the imputation we discarded geno- 
typed SNPs that were not present in the imputa- 
tion panel. For the GWAS, we added such SNPs 
back (if they passed quality control), for a total of 
7,428,049 SNPs (7,378,897 imputed and 49,152 geno- 
typed). To filter SNPs whose imputation results had 
changed over time, we performed an ANOVA test 
for frequency differences across batches. The quality 
control criteria for imputed SNPs were batch effects 
p- value at least 10~^^, average across batches of 
at least 0.5, and minimum across batches of at 



least 0.3. For genotyped SNPs, we required a MAF 
of at least 0.001, a Hardy- Weinburg p- value of at least 
10 and a call rate at least 0.9. 

Statistical analysis 

In order to minimize population substructure while 
maximizing statistical power, the study was limited 
to individuals with European ancestry. Ancestry 
was inferred from the genome-wide genotype data 
and principal component analysis was performed as 
in [44,47]. The cohort was filtered by relatedness to 
remove participants at a first cousin or closer relation- 
ship. More precisely, no two participants shared more 
than 700 cM of DNA identical by descent (IBD; ap- 
proximately the lower end of sharing between a pair 
of first cousins). IBD was calculated using the meth- 
ods described in [48]. The genomic control infiation 
factor was 1.156; all p- values are adjusted for this 
inflation factor. 

The GWAS was performed using likelihood ratio 
tests for the linear regression 

5 

i=l 

of carsickness on genotype, age, sex and 5 princi- 
pal components of genetic ancestry. Genotypes were 
coded as a dosage from 0-2 (counting the estimated 
number of minor alleles present for imputed SNPs) 
or as a count 0, 1, or 2 (also number of minor alleles, 
for genotyped SNPs). Significant SNPs were grouped 
into regions with at least 500kb between pairs of sig- 
nificant SNPs; the SNP with the lowest p- value in 
each region was chosen to be the index SNP. As some 
of the regions were under 1Mb apart, a joint regres- 
sion with all index SNPs was run to make sure that 
they all represented independent signals. 

Partial correlation between carsickness C and a 
phenotype Y were computed by computing the corre- 
lation between the residuals produced by regressing 
both C and Y on age, sex, and five principal compo- 
nents, using linear regression even if Y was a binary 
trait. We did not attempt to quantify the significance 
of these regressions nor any causality. 
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Tests of SNPs associated with motion sickness 
against other correlated traits were done using lo- 
gistic or linear regressions as appropriate with the 
same covariates as in the GWAS (except for morning 
sickness, which dropped sex). The phenotypes stud- 
ied (PONV, migraines, hay fever, altitude sickness, 
morning sickness, and vertigo) were all case control 
except for morning sickness, which was scored on a 5 
point scale: None; Mild (occasional bouts of queasi- 
ness or nausea, did not require treatment); Moder- 
ate (nausea and some vomiting, but did not require 
treatment); Severe (Severe nausea and vomiting that 
required treatment); Very severe (requiring hospital- 
ization and intravenous fluid (IV) therapy). 

Enrichment analysis using GREAT was conducted 
on all regions with an index SNP with p < 10~^, 
where regions were enforced to be 500kb apart. Win- 
dows of 500kb on eather side of each index SNP were 
uploaded into GREAT using default settings. 
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rsll713169 


-0.058 


-0.047 


1 


.233 


0. 


.415 


rsl782032 


-0.027 


-0.035 


0 


.774 


0. 


.434 


rs7170668 


0.038 


0.031 


1 


.224 


0. 


.482 


rs60464047 


-0.047 


-0.040 


1 


.182 


0. 


.601 


rs705145 


-0.049 


-0.053 


0 


.910 


0, 


.630 


rs2153535 


0.043 


0.048 


0 


.903 


0. 


.633 


rs9834560 


-0.043 


-0.039 


1 


.098 


0. 


.696 


rs2318131 


0.040 


0.036 


1 


.106 


0. 


.706 


rs7957589 


-0.044 


-0.049 


0 


.917 


0. 


.776 


rs4076764 


0.034 


0.032 


1 


.067 


0. 


.836 


rs4343996 


0.033 


0.034 


0 


.949 


0. 


.856 


rs2360806 


0.046 


0.048 


0 


.976 


0. 


.931 


rsl7515225 


0.032 


0.031 


1 


.018 


0. 


.954 



Table SI: Sex-specific effects for significant SNPs. Estimated effect sizes for females, males, and their 
ratio, p-values are for a non-zero interaction term between sex and genotype. SNPs are sorted by increasing 
p- value. 
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Figure S2: Region plots for genome-wide significant regions. SNPs are colored by with the index 
SNP (which is labeled). Circles are imputed SNPs, plus signs are genotyped SNPs. 
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